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ABSTRACT 

We study the impact of outflows driven by active galactic nuclei (AGN) on galaxy formation. Outflows 
move into the surrounding intergalactic medium (IGM) and heat it sufficiently to prevent it from con- 
densing onto galaxies. In the dense, high-redshift IGM, such feedback requires highly energetic outflows, 
driven by a large AGN. However, in the more tenuous low-redshift IGM, equivalently strong feedback 
can be achieved by less energetic winds (and thus smaller galaxies). Using a simple analytic model, we 
show that this leads to the anti-hierarchical quenching of star- formation in large galaxies, consistent with 
current observations. At redshifts prior to the formation of large AGN, galaxy formation is hierarchical 
and follows the growth of dark-matter halos. The transition between the two regimes lies at the z « 2 
peak of AGN activity. 

Subject headings: galaxies: evolution - quasars: general 



1. INTRODUCTION 

Downsizing is not a recent trend. For the past 10 billion 
years, the characteristic stellar mass of forming galaxies 
has been decreasing, as demonstrated by a variety of ob- 
servations. Space-based near-ultraviolet (NUV) measure- 
ments show that the typical star-formation rate in galax- 
ies was over an order of magnitude higher at z = 3 than 
at z = (Arnouts et al. 2005). Ground-based optical 
and near-infrared (NIR) searches indicate that the largest 
galaxies were already in place by z « 2, while smaller ones 
continued to form stars at much lower redshifts (Fontana 
et al. 2004; Glazebrook et al. 2004; van Dokkum 2004; 
Treu et al. 2005). And both optical and X-ray surveys 
detect a steady decrease in the characteristic luminosity 
of active galactic nuclei (AGN) below z « 3 (Pei 1995; 
Ueda et al. 2003) , which is likely to parallel the formation 
history of early- type galaxies {e.g. Granato et al. 2001). 

Yet despite such widespread evidence, galaxy downsiz- 
ing (Cowie et al. 1996) was unexpected. The A Cold Dark 
Matter (ACDM) model, while in spectacular agreement 
with observations {e.g. Spergel et al. 2003), is a hierarchi- 
cal theory, in which gravitationally-bound subunits merge 
and accrete mass to form ever-larger objects. Superim- 
posed on this DM distribution is the observed baryonic 
component. In canonical models, this gas falls into poten- 
tial wells, is shocked heated to the viral temperature, and 
must radiate this energy away before forming stars (Rees 
& Ostriker 1977; Silk 1977). The larger the structure, the 
higher its virial temperature, and the longer it takes to 
cool. Also, in most semi-analytical models, the dynami- 
cal time-scale is adopted for the rate of star formation via 
the Schmidt-Kennicutt law, and this time-scale increases 
systematically with halo mass. Thus larger galaxies were 
expected to form later not only because they occur in later- 
forming dark matter halos, but because the cooling and 
star formation times within such halos are longer. 

Two recent lines of inquiry have modifled this pic- 



ture. Theoretically, simulations have shown that cooling is 
greatly enhanced by gas inhomogeneities, allowing galax- 
ies to form with stellar masses over ten times larger than 
observed (Suginohara & Ostriker 1998; Dave et al. 2001). 
In fact, virializing shocks may even be completely absent 
in smaller halos (Birnboim & Dekel 2003). Observation- 
ally, it has been discovered that the main heating source 
for groups and small galaxy clusters is nongravitational 
{e.g. Arnaud & Evrard 1999). As these are the gaseous 
halos that were slightly too massive to form into galaxies, 
the clear implication is that nongravitational heating also 
played a key role in the history of large galaxies. 

Furthermore, the ~ 1 keV per gas particle necessary to 
preheat the intracluster medium (IGM) appears to exceed 
the energy available from supernovae {e.g. Valageas & Silk 
1999; Wu et al. 2000; Kravtsov & Yepes 2002). At the 
same time, active galactic nuclei (AGN) are observed to 
host high-velocity outflows with kinetic luminosities that 
may equal a signiflcant fraction of their bolometric lumi- 
nosity {e.g. Chartas 2002; Pounds et al. 2003; Morganti 
2005; see however Sun et al. 2005). Such outflows would 
naturally preheat the IGM to the necessary levels (Roy- 
chowdhury et al. 2004; Lapi et al. 2005), placing this gas 
at an entropy at which it can never cool within a Hub- 
ble time (Voit & Brian 2001; Oh & Benson 2003). This 
results in a feedback-regulated picture of galaxy forma- 
tion fundamentally different from the canonical approach 
(Scannapieco & Oh 2004; Binney 2004; Di Matteo et al. 
2005). 

In this letter, we explore the connection between this de- 
veloping picture and the widespread observations of down- 
sizing, structuring our investigation as follows: In §2 we 
develop a simple cooling model for galaxy formation, which 
illustrates how hierarchical formation naturally arises in 
the canonical approach. In §3 we modify our model to 
include AGN feedback, and show how this leads to a radi- 
cally different, anti-hierarchical history. We compare these 
results with recent observations and conclude with a short 
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discussion in §4. 

2. A COOLING-REGULATED MODEL OF GALAXY 
FORMATION 

Throughout this study, wc assume a ACDM model with 
parameters h = 0.65, ilo = 0.3, fl\ = 0.7, fib = 0.05, 
(T8 = 0.87, and n = 1, where h is the Hubble constant in 
units of 100 km Mpc~^, fio, f^A, and Qfc are the total 
matter, vacuum energy, and baryonic densities in units of 
the critical density, (t| is the variance of linear fluctuations 
on the 8/i^^Mpc scale, and n is the "tilt" of the primordial 
power spectrum [e.g. Spergel et al. 2003). The Eisenstein 
& Hu (1999) transfer function is adopted in all cases. 

We begin with a model of cooling-regulated galaxy for- 
mation and adopt as simple an approach as possible, to 
highlight the key physical issues. Unlike more sophisti- 
cated approaches {e.g. Kauffmann et al. 1993; Somerville 
& Primack 1999; Benson et al. 2000), we do not attempt to 
trace the detailed history of a sample of DM halos through 
the use of statistical "merger trees." Instead, we use the 
fit described in van den Bosch (2002), which gives the av- 
erage history of a DM halo with a present total mass of 
Mo as 
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where v = 1.34 -h 1.86 log(l -h-Zfit) -0.03 log(Mo/lOl2J^^0), 
M{z) is the mass of this system at some redshift z, and 
Zfit is defined implicitly through the relation _D(zfit)~"'^ ~ 
1 -t- 0.40v'o-2(0.254Mo) - ^^(Mo), with D being the Un- 
ear growth factor and (t^(M) being the variance of linear 
fluctuations within a sphere enclosing a total mass M. 

Following the canonical picture, we assume that at each 
redshift gas flows onto the growing halo along with the 
dark matter, is shocked to its virial density and temper- 
ature at that redshift [180 times the mean density and 
7.2 X lO'^ (Af(z)/10i2MQ)2/3 (1 + z) K, respectively], and 
is finally added to the pool of cold gas that forms stars 
after a delay of 

Uoi = 1.3 Te K)"' C-^K-_l:,{n) Myrs. (2) 

Here (rie) is the average number of electrons per cm'' in 
the shocked gas, C = (n^) (ne)~ is a "clumping factor," 
which accounts for inhomogeneities, Tg is the temperature 
in units of 10^ K, and A_23 is the radiative cooling rate 
of the gas in units of 10~^^ ergs cm^ s~^. Near 1 keV 
K. 10^ K, A_23 is roughly a constant but in general it is a 
function of temperature and metallicity. For simplicity we 
adopt the Sutherland & Dopita (1993) equilibrium values 
for this rate at a fixed metallicity of Q.IZq, the level of 
pre-enrichment of G dwarf stars in the solar neighborhood 
(Ostriker & Thuan 1975) and in large neighboring galax- 
ies (Thomas et al. 1999). Finally in keeping with previ- 
ous semianalytical models, we ignore gas inhomogeneities, 
taking C = 1. Adopting a larger value, consistent with 
simulations, would result in the same overall trends, but 
with galaxies that were much too large. 

We then compute the mass of cooled g cLS cLS cL func- 
tion of the total mass and redshift of observation, as il- 
lustrated in Figure 1. Here we see the classic behavior 
first described in Rees & Ostriker (1977) and Silk (1977). 



As their cooling times are short, Mcooi/Mtot ~ ^b/^o in 
small objects, while for larger masses with longer cooling 
times Mcooi/Aftot <S ^b/^o- In this picture the maximum 
.^cooi(2;) is the largest gas mass that has time to virial- 
ize and cool, a value that increases with time. Roughly 
one can estimate this scaling as fcooia = tf {1 + Zf )~^^'^ 
where tf and Zf are the final time and redshift at which the 
gas cools while ^cooi.i is computed at the initial redshift at 
which the gas falls on the halo. For the accretion histories 
described by eq. (1), Tvir depends only weakly on redshift, 
which allows us to rewrite this as Mgooi 2S (l + -z/)~^^'^Ivir- 
In Figure 1 we also compute the average redshift of for- 
mation of the resulting stars, adopting both a mass aver- 
age, that is zm = z{tM) where 



tM 



(3) 



and a U-band luminosity-weighted average, that is zlb = 
zitLs) where 



/o°"°dtTMeooiT^^(^obs-^0 



(4) 



and TB{t) is the i?-band mass to light ratio of a popu- 
lation of stars with an age t, assuming a Salpeter initial 
mass function as computed by Bruzual & Chariot (2003). 
In all cases zm and zlb are very close to the observed 
redshifts, meaning that in objects of all masses and red- 
shifts, most of the stellar mass and luminosity comes from 
recently-formed stars. 



3. 



A FEEDBACK-REGULATED MODEL OF GALAXY 

FORMATION 



We now modify our model to include AGN feedback, 
adopting an approach that parallels Scannapieco & Oh 
(2004), itself an extension of Wyithe & Loeb (2003). These 
papers showed that imposing an Mbh oc relationship 
betwecin black-hole mass and halo circiilar velocity (Fer- 
rarcsc 2002) and assuming that these objects shine at their 
Eddington luminosity for a fraction of the dynamical time 
after each major merger gives a good fit to the observed 
AGN luminosity function. In particular, each merger with 
a mass ratio less than 4:1 was associated with an AGN 
with a total bolometric energy (in units of 10^" ergs) of 



£;60,boi = 23Mf(^^„„i(l + z) ergs 



(5) 



where Mi2,cooi is the mass of cooled gas in miits of lO^^M©. 

Each AGN was then taken to host an outfiow, with 
a total kinetic energy equal to Ck « 0.05 of the total 
bolometric energy. Adopting a Sedov- Taylor model this 
resulted in a postshock temperature of T's(-R) = 8.8 x 
10^ek£^60,boiAfi2|gas(^)> where Mi2,gas is the total gas mass 
contained with a radius R in imits of 10^"^ Mq. Our simple 
model here does not distinguish between merger events, 
and so instead we relate the change in temperature at a 
radius R to the change in cooled mass by replacing Mi2,cooi 
with dMi2,cooi : 



dTsjR) 
dM 



12, cool 
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Apart from this additional contribution to the tempera- 
ture, our model exactly parallels that in §2, with two minor 
modifications. First, if a parcel of gas remains at a temper- 
ature above the final virial temperature of the halo for a 
time longer than its dynamical time, we assume it escapes 
from the gravitational potential. Secondly, we choose the 
clumping factor to be consistent with simulations, which 
show that roughly 40% of the gas cools onto galaxies by 
2; = in the absence of feedback (Dave 2001; Balogh et al. 
2001). This requires us to raise Mcooi (without feedback) 
by a factor ^ 30, which corresponds to C « 3. 

With these modifications, and adopting a fiducial value 
of eic = 0.05, we constructed plots of Mcooi, zm and zlb 
as in §2. These are shown in the lower panels of Figure 
1. As in our cooling model, Mcooi /Aftot ~ ^b/^a, in small 
halos, where feedback is weak. At large masses, however, 
feedback causes Mcooi to evolve radically differently. Now 
the maximum scale of star-forming galaxies drops at lower 
redshifts. This is the direct result of the decrease in the 
cosmological gas dcinsity and the fact that at lower red- 
shifts smaller AGNs (and hence smaller galaxies) can heat 
this gas to temperatures that will not cool within a dy- 
namical time. The scaling of this "quenching threshhold" 
(Faber et al. 2005) can be estimated by setting the dy- 
namical time (oc n^^/^) equal to the post-AGN feedback 
cooling time given by eqs. (2) & (5). For the « 1 keV 
objects we are interested in this gives Mcooi SS (1 + -z)"^^''- 

Galaxies with the largest Mcooi values (f« 3 x IO^^Mq) 
are already in place at 2: = 3 and continue to fall into 
larger halos at lower z without increasing in stellar mass. 
Somewhat smaller Afcooi ~ lO^^Af© galaxies can recover 
from an AGN outflow at z = 1.5, but not at the lowest 
redshifts. Finally galaxies with Mcooi < IO^^Mq have can 
hold onto their gas, even at 2 = 0. Note, however, that our 
simple models do not include supernova feedback, which 
is likely to be important in galaxies with Mtot <, IQ^^Mq 
(Martin 1999). This mechanism should behave qualita- 
tively differently, as it depends primarily on momentum 
transfer rather than heating (Thacker et al. 2002). Note 
also that we do not account for the possibility of initial 
positive AGN feedback before quenching occurs, as may 
be required by observations of efficient star formation in 
massive high-redshift galaxies (Silk 2005). 

The anti- hierarchical trend in our model is also apparent 
from the mass and i?-band luminosity-weighted ages. In 
small galaxies, most of the stars are relatively young and 
were formed at redshifts close to the redshifts of observa- 
tion. Stars in large galaxies, however, are much older, and 
their ages increase as a function of Mcooi- 

4. COMPARISON WITH OBSERVATIONS AND DISCUSSION 

The features of the simple models developed in §2 and §3 

can be directly compared with observations. As a measure 
of the typical star formation rate, we compute the absolute 
AB magnitude at 1500 A for our model galaxies, assuming 
that 50% of the cooled gas is in the form of stars, adopting 
the Bruzual fc Chariot (2003) population synthesis mod- 
els, and ignoring dust. We do not attempt to construct 
a luminosity function of galaxies, since doing so would 
require us to use a full merger-tree formalism. Rather, 
for both cooling and AGN feedback models, we impose a 
halo mass limit of :^(Mnl, z) = 5c a{M)-'^D{z)''^ < 2.5 at 



each redshift, to exclude excessively rare peaks. For the 
CDM power spectrum at these redshifts, A/nl cx; {1 + z)~^, 
such that the corresponding virial temperature Tvir,NL 

il + z)-'^A4i\ 

We then find the most luminous galaxy at 1500 A within 
the remaining halos, and plot its absolute magnitude as a 
function of redshift in the upper panel of Figure 2. This is 
compared with the observed evolution of Afj*5QQ to z ^ 3, 
which is divided into two regions. At high redshift, it 
brightens along with the nonlinear mass scale over the dy- 
namical time. At low redshift, the most luminous galaxies 
no longer lie in the largest halos. Rather Mfgoo fades along 
with the AGN quenching mass scale over tdyn- The redshift 
at which these two scales cross marks a distinct transition 
between hierarchical and anti-hierarchical growth, which 
occurs at the peak of AGN activity {e.g. Pei 1995; Ueda 
et al. 2003). 

The shaded regions in this figure show the impact of 
varying Ck, which shifts the peak value of M^^qq. Modi- 
fying our other assumptions (Mnl, cold fraction in stars, 
etc.) would result in similar changes. In particular, dust 
would cause M^c^qq to saturate at z « 2 due to the larger 
attenuation factors of large starbursts {e.g. Wang & Heck- 
man 1996; Adelberger & Steidel 2000; Martin et al. 2005). 
Our point is not to fit these particular parameters to the 
data, howevever, but rather that substantial high-redshift 
brightening and lower-redshift fading of M^^qq are general 
and unavoidable features of AGN-regulated galaxy forma- 
tion. 

Coincidentally, this sharp peak in Mj^gQ,, is now starting 
to be seen observationally. Whereas it has been known for 
several years that M* brightens monotonically from z ^ 
to 2; f« 2 in the UV and bluer optical bands {e.g. Gabasch 
et al. 2004; Arnouts et al. 2005), recent measurements at 
z « 5 - 6 (Ouchi et al. 2004; Bouwens et al. 2005) now 
also show a downturn relative to lower redshifts, implying 
that MfgQo is most luminous around 2 « 3. Remarkably, 
this is very close to the peak in AGN activity, suggesting 
that z « 2 — 4 represents the key epoch for gas accretion 
onto massive systems. 

Turning our attention to the evolution of the maximum 
stellar mass, as shown in the bottom panels, similar trends 
are apparent. While the maximum mass in the cooling 
model increases monotonically, the AGN feedback model 
is divided into two regions: one at high redshift in which 
the maximum stellar mass increases along with the non- 
linear mass scale, and one at low redshift in which the 
maximum stellar mass stays fixed as the scale at which 
new galaxies are forming becomes smaller. The general 
properties of cooling and AGN feedback models are sum- 
marized in Table 1. 

Finally, we point out that not only the temporal, but the 
spatial distribution of AGN in our feedback model is sug- 
gestive of recent observations. Using over 20,000 galaxies 
from the 2dF QSO Redshift Survey, Groom et al. (2005) 
have studied the spatial clustering of QSOs near the char- 
acteristic scale in the optical luminosity function. They 
measure bias values which, when converted into i^(z) val- 
ues using standard expressions (Mo & White 1996), evolve 
from 2.5 ± 0.2 at z = 2.48 down to 1.1 ± 0.2 at z = 0.56. 
The characteristic scale of low-redshift AGN is downsizing 
from the rarest to the most common objects, spreading 
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Table 1 

Properties of Galaxy Formation Models 



Model 


Characteristic Characteristic Characteristic z = Trends 
A^cooi at high- 2; Mcooi at low-^ SFR at low- 2; 


Cooling 


Tracks Mnl Cooling-regulated Mcooi/tdyn Bigger galaxies are 

°S (1 + z)~^ RS (1 + 2;)"''^^ roughly constant young and blue 


AGN 
Fccnlback 


Tracks ilfxL F(X!dbac;k-rc^guIat{xl i\fcooi/'dvn Biggc^r galaxi(^s are 

+ ;^ (1 + ,:)■■''' 2S(1 + ■-)''' ^ okltT and rockier 



the heated gas that extinguishes the formation of galaxies 
to this day. 
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Fig. 1. — Mass of cooled gas in models regulated by cooling (top panels) and AGN feedback (bottom panels). In both cases the left panels 
compare the cooled gas mass with the total mass, the center panels give the mass-averaged star-formation redshift, and the right panels give 
the luminosity-weighted (B-band) average star-formation redshift. In all panels the solid, long-dashed, and short-dashed lines correspond 
to galaxies observed at redshifts of 0, 1.5, and 3 respectively. In the cooling-regulated models (§2), larger galaxies continue to form at late 
times, and all galaxies are dominated by young stars. In the models regulated by AGN feedback (§3), the maximum scale of forming galaxies 
decreases with time, and large galaxies are accreted into bigger halos without forming new stars. When observed at low redshift, the stars in 
these objects are older. 
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Fig. 2. — Top: Evolution of JV/J'jqq, the characteristic scale for the instantaneous star formation rate in galaxies. Lines are from our cooling 
model (dashed) and ej; = 0.05 AGN feedback model (solid) compared with measurements by Arnouts et al. (2005: solid points) and Steidel 
et al. (1999: open square). Bottom: Evolution of the characteristic stellar mass as compared to measurements by Fontana et al. (2004). 
Lines are as in the upper panel. In both panels, the shaded regions correspond to AGN feedback models with ej^ varied from 0.025 to 0.1. 
The AGN feedback model is divided into two regimes: a high-redshift ( "coalescence" ) regime in which the characteristic scale of star-forming 
galaxies increases with the nonlinear mass scale, and a low-redshift ("feedback") regime in which this characteristic scale decreases as a result 
of quenching from ever more efficient feedback. 



